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基于 二 分 类 Logistic 回归 模型 的 太行 山 丘 陵 区 
县 域 耕 地 资源 潜力 估算 - 
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摘 要 : 耕地 红线 划 定 与 人 -地 资源 矛盾 日 益 突 出 背景 下 ， 耕 地 资源 潜力 的 研究 与 开发 日 显 重要 。 ee 
积 近 2/3 分 布 在 山区 ， 因此 山区 耕地 资源 的 合理 开发 利用 及 其 资源 潜力 的 研究 尤为 重要 。 本文 以 华北 地 区 的 太 
行 山 为 研究 区 域 , 选择 耕地 占 比 和 资源 潜力 最 大 的 丘陵 区 典型 县 一 一 河北 省 
响 耕 地 资源 潜力 的 基本 生态 要 素 , 包括 5 个 地 形 要 素 和 8 个 直接 气象 要 素 或 由 气象 要 素 计算 得 到 的 间接 气象 
要 素 ， 引 入 二 分 类 Logistic 回归 分 析 方 法 ,运用 偏 最 大 似 然 估计 向 前 引入 法 的 拟 合 方法 ,筛选 提取 影响 耕地 资 
源 潜力 的 关键 生态 要 素 ; 由 模型 参数 Waldx2 统计 量 分 析 影 响 耕 地 资源 潜力 的 关键 生态 要 素 的 贡献 率 排序 ; 由 
模型 参数 回归 系数 6 分 析 耕 地 资源 潜力 与 生态 要 素 的 相关 关系 ; 由 模型 参数 发 生 比率 OR 分 析 量 化 关键 生态 
要 素 对 耕地 资源 潜力 的 影响 ， 最 终 建立 Logistic 回归 模型 . 基于 此 模型 ,在 GIS 软件 中 得 到 井 取 县 耕地 资源 潜 
力 分 布 图 , 进而 估算 出 县 域 耕 地 资源 潜力 .研究 结果 表明 : 13 个 影响 井 取 县 耕地 资源 潜力 的 基本 生态 要 素 中 8 个 
为 关键 生态 要 素 ; 关键 生态 要 素 中 地 形 要 素 配 置 比 气象 要 素 配 置 更 为 重要 ; 年 平均 气温 和 寒冷 指数 与 耕地 资源 
潜力 呈 负 相关 关系 ,其 余生 态 要 素 则 呈正 相关 关系 ; 由 回归 模型 估算 出 井 取 县 具备 垦殖 为 耕地 资源 的 土地 面积 
为 60 400 hm”, 而 根据 晃 感 影像 解 译 结果 得 出 的 现 有 耕地 资源 为 45 600 hm”， 由 此 井 取 县 尚 共 有 14 800 hm 的 后 
备耕 地 资源 , 相当 于 现 有 耕地 面积 的 32.5%, 这 说 明 在 不 考虑 垦殖 所 带 来 的 可 能 负 效 应 的 前 提 下 ， 井 险 县 具有 和 较 
大 的 后 备耕 地 资源 开发 潜力 ， 该 结论 为 井 降 县 后 备耕 地 资源 的 开发 与 可 持续 利用 提供 了 理论 依据 。 

关键 词 : 太行 山 丘 陵 区 ; 二 分 类 Logistic 回归 模型 ; 生态 要 素 ; 耕地 资源 潜力 ; 后 备耕 地 资源 
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Abstract: Research on the potential of cultivated land resources in mountainous area is particularly important for the 
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sustainable development and utilization of arable lands. Using Jingxing County of Hebei Province, a typical hilly county in 
Taihang Mountain, as a case Study, the potential for reserved cultivated land in mountain regions was estimated. A total of 13 
ecological factors influencing the potential of cultivated land resources in hilly areas were used in the analysis 一 3 terrain 
factors and 8 climate or climate-related factors. A binary Logistic regression model was built on key ecological factors to 
analyze the potential of arable land resources in mountain regions in the case-study area. The distribution map of potential 
arable land resources in Jingxing County was draw in GIS environment based on the result binary Logistic regress and the 
potential arable land resources in the region analyzed from the map. The key ecological factors among 13 factors that 
influenced arable land resources were extracted using the same binary Logistic regression model. The order of the contribution 
of key ecological factors was determined using the model parameter (Waldx). The regression coefficient (8) was used to 
analyze the correlation between the key ecological factors and potential arable land resources in the region. Odds ratio (OR) 
was use to show correlation between the changes in configuration of key ecological factors and the changes in potential arable 
land resources. The results showed that there were 8 key ecological factors influencing cultivated land resources in the 
case-study area, with the importance order of relief > elevation > slope position index > trasp > coldness index > annual 
average temperature > drought index > slope. Terrain factors were more important indicators than the climate-related factors in 
estimation of cultivated land resource potential. Among 8 key ecological factors, annual average temperature and coldness 
index were both negatively and positively correlated with potential arable land resources, while other 6 factors were all 
negatively correlated with potential arable land resources. OR calculation showed that except drought index, the one unit 
change of the key ecological factors caused 1 time change in cultivated land resource potential. About 60 400 hm” of land 
resources Were available for cultivation in the Jingxing County. The analysis of SPOTS image showed 45 600 hm’ of existing 


cultivated land. Therefore, there were 14 800 hm’ 


of cultivatable land reserves in the county, accounted for 32.5% of the 
existing cultivatable lands in the region. This suggested that there still were a lot of cultivatable land reserves in Jingxing 
County. The results of this study provided theory basis for the development and utilization of cultivated land reserves in 
Jingxing County. 

Keywords: Hilly area of Taihang Mountain; Binary Logistic regression model; Ecological factor; Cultivated land resource 


potential; Cultivatable land reserve 


山地 约 占 全 球 陆地 总 面积 的 20%05， 世 界 约 
12% 的 人 口 居住 在 山区 户 。 我 国 是 一 个 多 山 国家 , 山 
地 面积 约 占 陆地 总 面积 的 70%，40% 多 的 人 口 生活 
在 山区 5 1。 耕地 是 人 类 赖 以 生存 的 根本 ， 我 国 耕 地 
中 有 近 2/3 分 布 在 山区 、 丘陵 和 高 原 ， 且 耕地 面积 在 
持续 减少 5l。 进 入 21 世纪 以 来 , 由 于 山区 退耕 还 林 
还 草 政 策 实施 ,山区 耕地 面积 逐渐 减少 ,而 社会 
经 济 持续 快速 增长 也 不 可 避免 的 要 占用 耕地 。 因 此 ， 
为 确保 “18 亿 亩 耕地 的 红线 ”， 研 究 利 用 后 备耕 地 资 
源 对 保障 耕地 动态 平衡 意义 重大 。 我 国 已 经 进行 了 
几 次 大 规模 的 后 备耕 地 资源 调查 ,数据 显示 ,后 备 
耕地 资源 中 面积 最 大 的 是 荒草 地 ， 占 总 后 备耕 地 资 
源 的 88.44%0]。 

当前 对 山区 耕地 资源 方面 的 相关 研究 主要 集中 
于 3 个 方面 : 一 是 对 山区 耕地 资源 时 空 分 布 特征 的 
研究 ,在 时 间 上 利用 遥感 影像 解 译 、 结 合 地 理 信息 
系统 (GIS) 空 间 分 析 功 能 研究 不 同 区 域 尺度 的 土地 
利用 (包括 耕地 ) 时 空 变 化 , 如 Ullah 等 外 研究 了 兴 都 
库 什 山脉 耕地 的 时 空 变 化 特征 ; 在 空间 上 耕地 资源 
的 分 布 格局 特征 主要 从 地 形 因 素 ( 高 程 、 坡度) 方面 进 
行 分 析 ; 谭 莉 梅 等 站 对 河北 省 太行 山区 域 耕 地 资源 
的 空间 分 布 从 海拔 、 坡 度 及 地 形 类 型 3 个 方面 进行 


了 分 析 研 究 。 二 是 对 山区 耕地 资源 的 适宜 性 评价 研 
究 , 胡 学 东 等 "以 自然 因素 为 主 的 8 个 评价 指标 评 
价 鄂 西北 山地 区 耕地 整治 的 适宜 性 ; Liu 等 "从 气 
候 、 水 文 、 地 形 、 土 壤 和 植被 5 个 方面 对 土地 适宜 
性 进行 评价 ,得 出 部 分 秦岭 山区 耕地 需要 退耕 。 三 
是 对 山区 后 备耕 地 资源 的 适宜 性 研究 ， 张 甘霖 等 1 
和 张 迪 等 上 分 别 评价 了 后 备耕 地 资源 的 自然 质量 适 
宜 性 和 经 济 供给 适宜 性 , 研究 中 所 用 的 后 备耕 地 资 
源 只 是 根据 土地 利用 分 类 数据 依据 耕地 后 备 资源 定 
义 中 提取 得 到 的 , 而 且 其 研究 对 象 也 不 仅 是 山区 。 
前 人 的 研究 中 ， 或 是 针对 现 有 耕地 资源 存量 的 研究 ; 
或 是 土地 利用 分 类 数据 、 偏 重 于 自然 质量 和 经 济 供 
给 适宜 性 ， 过 于 侧重 从 土地 管理 部 门 的 需求 角度 出 
发 。 因 此 ， 本 研究 重点 关注 耕地 资源 潜力 和 后 备耕 
地 资源 的 相关 性 ,在 考虑 现 有 耕地 存量 条 件 下 ， 综 
合 考 虑 山区 影响 耕地 资源 潜力 的 地 形 与 气象 等 生态 
要 素 , 通过 建立 模型 ， 精 准 估算 土地 可 垦殖 为 耕地 
资源 的 综合 潜力 及 其 分 布 , 并 与 现 有 耕地 存量 及 分 
布 比较 ， 从 自然 生态 条 件 和 战略 储备 上 分 析 探讨 山 
区 尚 具备 的 后 备耕 地 资源 开发 潜力 ， 这 个 方面 的 研 
究 目 前 尚 少见 报道 。 

张 迪 等 上" 指出 后 备耕 地 资源 中 , 面积 最 大 的 是 
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草地 ， 这 部 分 荒草 地 大 多 分 布 在 山地 丘陵 区 "1。 张 
贵 军 等 nl 按 不 同 地 貌 类 型 将 河北 省 划分 为 8 个 类 型 
区 ， 其 中 太行 山 丘 陵 区 耕地 后 备 资源 潜力 较 大 。 
此 ， 本 文 以 太行 山 为 研究 区 域 ， 选 择 耕 地 占 比 和 资 
源 潜力 最 大 的 丘陵 区 为 研究 对 象 ， 以 河北 省 井 降 县 
为 典型 县 域 研究 案例 ， 建 立 丘 陵 区 耕地 资源 潜力 与 
关键 生态 要 素 的 二 分 类 Logistic 回归 模型 。 基 于 此 
模型 ,在 GIS 软件 中 得 到 太行 山 丘 陵 区 井 降 县 耕地 
资源 潜力 分 布 图 ,进而 估算 出 县 域 耕 地 资源 潜力 ， 
并 与 根据 遥感 解 译 所 得 的 井 降 县 耕地 资源 信息 对 比 
分 析 , 为 井 取 县 后 备耕 地 资源 的 开发 、 政 府 决 策 和 
耕地 资源 可 持续 利用 提供 理论 依据 。 


1 研究 区 域 概况 与 研究 方法 
1.1 ”研究 区 概况 

井 取 县 位 于 河北 省 西部 ， 太 行 山东 苞 ，37°42'~ 
38°13"N, 113"48'~114?18'E， 是 100% 山 区 县 。 井 取 县 
总 面积 138 320 hm”， 其 中 耕地 22 867 hm ， 占 总 面 
积 16.5%。 全 县 总 人 口 32 万 。 井 降 县 属于 温带 亚 湿 
润 大 陆 性 气候 ， 海 拔 从 145~1 188 m， 降 水 量 456~ 
521 mm， 年 平均 温度 6~12 'C， 宇 10 CC 积温 3 273~ 
4 314 'C。 山 区 坡度 0*~58°。 井 取 县 在 整个 太行 山区 
的 位 置 如 图 1 所 示 。 
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1 井 取 县 在 太行 山区 位 置 图 
Fig.1 Map of Taihang mountains area and the location of the 


study area of Jingxing County 


1.2 ”数据 来 源 及 其 处 理 

为 便于 运算 ,所 有 的 数据 都 采用 GIS 的 栅 格 数 
据 类 型 ( 表 1), 并 且 在 90 mx90 m 的 分 辩 率 上 进行 
重 采样 ， 都 配 准 到 WGS84 坐标 系统 。 


表 1 遥感 影像 解 译 耕 地 资源 误差 矩阵 表 


Table 1 Error matrix of cultivated land resources of the 

remote-sensing image interpretation Grids 
遥感 影像 解 译 数据 类 型 
实际 数据 类 型 Data type of image interpretation 
Real data type 耕地 非 耕 地 总 和 
Cultivated land Non data type Total 
耕地 

1 8 14 2 036 

Cultivated land 
非 耕地 480 3 520 4000 

Non cultivated land 
总 和 Total 2 373 3 663 6 036 
1.2.1 ” 井 取 县 耕地 资源 


为 得 到 更 加 详细 、 准 确 的 耕地 资源 数据 ， 本 研 
究 利 用 2006 年 10 月 的 分 辩 率 为 2.5 m 的 SPOT5 遥 
感 数据 解 译 研究 现 有 耕地 。 解 译 中 将 土地 利用 分 为 
耕地 和 非 耕 地 两 种 类 型 (图 2)。 解 译 结果 精度 评价 采 
用 典型 区 域 野外 抽样 调查 ,在 解 译 为 耕地 的 区 域 随机 
选取 了 20 个 样 区 , 总 面积 为 12 725 m”, 包含 2036 个 
栅 格 单元 ， 耕 地 解 译 精度 为 93%; 在 解 译 为 非 耕 地 
的 区 域 随机 选取 40 个 样 区 ,总 面积 为 25 000 m2,， 包 
含 4000 个 栅 格 单元 , 非 耕 地 解 译 精度 为 88%， 总 的 
解 译 精度 为 90%( 表 1)。 


[下 ] 关 电 
| Arable land 


非 耕地 


Non-arable land 


2 2006 年 井 降 县 耕地 资源 分 布 图 
Fig.2 Distribution map of cultivated land of Jingxing 
County in 2006 
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1.2.2 ”影响 耕地 资源 潜力 的 生态 要 素 reos (jes ea-a0)| 
影响 耕地 资源 分 布 的 生态 因素 很 多 ， 根 据 胡 学 180) ? 
i 如 下 : Trasp ，Trasp 的 
东 等 "在 鄂 西 北山 地 区 耕地 整治 适宜 性 评价 中 选取 2 
的 影响 耕地 整治 的 因素 以 及 张 甘霖 等 (2 研究 的 中 国 。”” 值 域 为 [0, 1], 0 表示 接受 的 太阳 能 最 少 , 1 表示 接受 


后 备耕 地 资源 自然 适宜 性 评价 中 的 评价 因子 , 同时 
考虑 数据 的 可 获取 性 ， 本 研究 选取 了 5 个 地 形 要 素 ， 
包括 高 程 、 坡 度 、 转 换 坡 向 、 地 形 起 伏 度 和 坡 位 指 
数 , 8 个 气象 要 素 及 由 气象 要 素 计算 所 得 的 相关 要 素 ， 
包括 即 年 平均 气温 、 年 平均 降水 、 年 活动 积温 、 年 
生物 积温 、 温 暖 指数 、 寒 冷 指 数 、 湿 润 指 数 和 潜在 
蒸发 率 ( 表 2)。 


表 2 耕地 资源 及 影响 耕地 资源 潜力 的 生态 要 素 属 性 表 
Table 2 Attributes of cultivated land and ecological factors 
influencing cultivated land resources potential 


变量 数据 层 栅 格 类 型 单位 或 描述 
Variable Data Grid type Unit/describe 
开明 = 
耕地 二 分 类 
Sh non Arable land Binary 
variable 
恋 旦 -车 张 型 
自 变量 高 程 DEM 连 续 型 本 
(地 形 要 素 ) Continuous 
Independent 坡度 连续 型 
variable Slope direction Continuous 
(terrain factors) 照射 坡 向 连续 型 无 量 纲 
Trasp Continuous Non-dimensional 
地 形 起 伏 度 连续 型 本 
Relief Continuous 
坡 位 指数 a ， 
Topographic Continuous 
position index 
自 变 量 年 平均 温度 连续 型 和 
(气象 要 素 ) Annual average ti 
Dependent temperature 
| variable 年 平均 降水 连续 型 
(climate factors) Annual average g mm 
ee Continuous 
precipitation 
年 活动 积温 。 连续 型 | 
Active accumulate : C 
Continuous 
temperature 
年 生物 积温 
Biological 连续 型 CC 
cumulative Continuous 
temperature 
温度 指数 连续 型 C 
Warmth index Continuous 
寒冷 指数 连续 型 C 
Coldness index Continuous 
湿润 指数 连续 型 | 
a 。 mm C 
Humidity index Continuous 
干燥 度 连续 型 无 量 纲 
Drought index Continuous Non-dimensional 


高 程 (DEM) 数 据 来 自分 辩 率 为 90 m 的 美国 
SRTM 数据 (http://www.gscloud.cn )( 图 3A); 以 DEM 
为 原始 数据 ， 坡 度 、 坡 向 分 别 利用 ARCGIS10.2 中 的 
slopevaspect 命 令 直 接 得 到 各 个 栅 格 的 坡度 值 (图 3B) 
和 坡 向 值 。 照 射 坡 向 (Trasp) 能 够 直观 反映 接受 太阳 
能 多 少 ， 需 要 对 坡 向 (aspecb 进 行 一 定 的 转换 ， 公 式 


的 太阳 能 最 多 (图 3C)。 地 形 起 伏 度 指 地 面 某 一 确定 
距离 的 范围 内 最 高 点 与 最 低 点 的 高 程 差 ， 本 文选 择 
3x3( 单 位 : 栅 格 ) 的 矩形 区 域 基于 GIS 得 到 各 个 栅 格 
的 地 形 起 伏 度 (图 3D)。 坡 位 指数 是 用 来 计算 某 一 单 
元 在 整体 景观 中 的 位 置 ， 通过 计算 栅 格 单元 与 周 
围 栅 格 单元 的 平均 海拔 差 得 到 的 ， 本 文选 择 尺度 为 
20x20( 单 位 : 栅 格 ) 的 矩形 区 域 利用 GIS 软件 得 到 各 
个 栅 格 的 坡 位 指数 (图 3E)。 

气象 数据 是 通过 1961 一 2005 年 气象 站 的 观测 
得 到 的 ， 需 要 对 多 年 数据 平均 值 进行 空间 差 值 ， 从 
而 得 到 每 个 栅 格 单元 的 气象 要 素数 据 。 本 文 利用 全 
720 个 气象 台 站 的 气温 和 降水 数据 结合 高 程 、 经 
度 和 纬度 ， 建 立 气 温和 降水 的 多 元 回归 方程 ， 然 后 
利用 内 插 的 方法 得 到 每 个 栅 格 的 日 气温 和 降水 数 
据 。 年 平均 气温 和 年 平均 降水 利用 每 个 栅 格 的 日 气 
温 、 日 降水 数据 计算 得 到 (图 3F-G)。 年 宇 10 人 活动 
积温 通过 ARCGIS10.2 中 的 raster calculator 命令 计 
算得 到 (图 3H)。 年 生物 积温 是 指 >0 ‘CC、<30 的 月 
平均 气温 (当月 平均 气温 >30 ‘CC 时, 取 值 为 30 'C) 之 
和 的 平均 值 (图 3D。 温 暖 指 数 是 指 一 年 中 月 平均 气温 
下 5 'C 以 上 的 月 份 的 月 平均 气温 与 5 它 只 差 的 累计 值 
(图 3j)。 塞 冷 指数 是 指 月 平均 气温 小 于 5 它 的 月 份 的 
月 平均 气温 与 5 只 差 的 累计 值 (图 3K)。 湿 润 指数 
是 降水 量 与 温暖 指数 的 比值 (图 3L)。 干燥 度 是 年 洪 
在 蒸发 量 与 降雨 量 的 比值 ， 本文 的 年 潜在 蒸发 量 根 
据 Holdridge 经 验 公 式 计 算得 到 ， 即 年 潜在 蒸发 量 等 
于 58.93 乘 以 年 生物 积温 (图 3M)。 以 上 气象 相关 要 
素数 据 都 是 根据 日 气温 和 日 降水 数据 利用 GIS 的 栅 
格 运 算得 到 的 。 
1.3 ”二 分 类 Logistic 回归 模型 
1.3.1 ”模型 原理 

二 分 类 Logistic 回归 分 析 是 针对 因 变 量 为 二 分 
类 的 非 线 性 回归 统计 方法 。 本 研究 中 耕地 ( 力 是 二 分 
类 因 变 量 。 当 土地 利用 为 耕地 时 芒 1， 否则 ,六 0。 
假设 影响 土地 利用 是 否 能 成 为 耕地 的 是 13 生态 要 
素 ， 即 自 变量 C 思 (=1，2，…, 及， 土地 利用 为 耕地 
(天 1) 的 概率 为 P。 二 分 类 Logistic 回归 模型 如 式 
(0D2 29: 

EXP(o+AXI+AAXa + + Pr Xr) 


1 
1+EXP( + BX1+PX; + + PrXi) 
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roger (7 -TjnfP/0-A)]- 四 


Xi 


Pot+BX1t+PB Xs + + PX 
式 中 : p;(i=0, 1, 2,…, 有 为 回归 系数 。P/(1-P) 为 事件 
的 发 生 比 率 (odds ratio)。 
1.3.2 ”模型 验证 
本 研究 用 HL 指标 检验 二 分 类 Logistic 回归 模型 
拟 合 优 度 。 当 HL 指标 统计 显著 表示 模型 拟 合 不 好 。 相 
反 , 当 HL 指标 统计 不 显著 表示 模型 拟 合 好 。HL 指标 
是 一 种 类 似 于 皮尔 逊 ?统计 量 的 指标 , 其 公式 如 下 5 
& (LSND) 
WN RP) , 
式 中 : G 为 分 组 数 ，G 科 10; ,为 第 g 组 中 的 案例 数 ; 
了】 为 第 g 组 事件 的 观测 数量 ; Ps 为 第 g 组 预测 事件 概 
率 ; NvP. 为 预测 数 。 
1.4 二 分 类 Logistic 回归 模型 建立 
1.4.1 样本 选取 
样本 选取 通过 Arcgis 中 的 随机 取样 命令 来 实现 ， 
首先 将 井 取 县 土地 利用 分 为 耕地 和 非 耕 地 两 个 栅 格 
图 ,利用 ArcgisWorkstation 中 的 rand 命令 对 所 有 的 
井 取 县 进行 随机 编码 。 然后 耕地 取 编 码 数值 小 于 0.4 
的 栅 格 ， 非 耕地 选择 编码 数值 小 于 等 于 0.2 的 栅 格 
作为 掩 膜 (也 就 是 选取 了 井 取 县 耕地 面积 的 40% 和 
井 取 县 非 耕 地 面积 的 20% 作 为 样本 )， 共 选取 的 样 点 
数 为 46 083 个 ， 其 中 耕地 的 样 点 数 为 22 254 个 ,， 非 
耕地 的 样 点 数 为 23 829 个 , 不 相等 的 抽样 比例 不 会 
影响 解释 变量 在 Logistic 回归 模型 中 的 系数 估计 ， 
但 是 会 影响 模型 的 常数 项 7， 为 保证 耕地 与 非 耕 
地 样 点 数 的 基本 一 致 ， 对 耕地 和 非 耕 地 分 开 取样 . 


选择 不 同 的 编码 数值 界线 (因为 山区 耕地 面积 较 小 )。 
最 后 , 应 用 GIS 软件 ,基于 井 取 县 关键 生态 要 素 栅 
格 图 提取 井 取 县 内 耕地 和 非 耕 地 样 点 所 对 应 的 13 
个 关键 生态 要 素 的 值 。 
1.4.2 ”模型 建立 

将 因 变 量 耕 地 资源 (耕地 / 非 耕 地 ) 和 相对 应 的 自 
变量 (13 个 关键 生态 要 素 值 ) 输 入 到 SPSS 软 件 中 , 利 
用 软件 中 二 分 类 Logistic 回 归 模 型 ， 选 择 偏 最 大 似 然 
估计 向 前 引入 法 的 拟 合 方 法 将 P 值 小 于 0.05 的 变量 
引入 回归 过 程 ,这 种 方法 中 ,被 选 自 变量 依次 进入 
回归 模型 ， 首 先 引 入 与 因 变 量 最 大 偏 相 关 且 符合 引 
入 标准 的 自 变量 ,在 引入 第 1 个 变量 后 ， 再 引入 第 2 
个 与 因 变 量 最 大 偏 相 关 并 且 符 合 引 入 标准 的 自 变量 ， 
依次 类 推 。 当 无 自 变 量 符合 引入 标准 时 ， 中 止 回归 
过 程 。 采 用 此 方法 ， 从 13 个 生态 要 素 因 子 中 选择 出 
影响 山区 耕地 资源 的 关键 生态 要 素 8 个 ， 引入 的 顺 
序 依次 为 地 形 起 伏 度 、 高 程 、 坡 位 指数 、 转 换 坡 向 、 
寒冷 指数 、 年 平均 气温 、 干 燥 度 、 坡 度 。 
1.4.3 ”模型 验证 

由 表 3 可 知 ,模型 的 HL 指标 为 2.532， 其 概率 为 
0.98， 统 计 不 显著 ,说 明 耕 地 资源 可 能 性 分 布 的 二 
分 类 Logistic 回 归 模 型 很 好 地 拟 合 了 数据 。 
2 结果 与 分 析 
2.1 山区 耕地 资源 潜力 的 关键 生态 要 素 

偏 最 大 似 然 估 计 向 前 引入 法 的 拟 合 方法 ， 从 13 
个 生态 要 素 因 子 中 选择 了 影响 山区 耕地 资源 的 关键 
生态 要 素 8 个 ,引入 的 顺序 依次 为 地 形 起 伏 度 、 高 
程 、 坡 位 指数 、 转 换 坡 向 、 寒 冷 指 数 、 年 平均 气温 、 
干燥 度 、 坡 度 。 由 模型 参数 Waldx 统计 量 ( 表 3) 可 知 ， 


表 3 耕地 资源 ( 因 变 量 ) 及 其 关键 生态 要 素 ( 自 变量 ) 的 二 分 类 Logistic 回归 模型 参数 及 显著 水 平 
Table 3 Parameters and significance ofbinary Logistic regression model of cultivated land (dependent variable) and its key 
ecological factors (independent variables) 


一 | 2 一 \ 作 、 口 区 人 7 三 旦 荔 ， 水 平 > < 
自 变 量 回归 系数 (8) 标准 误差 (S.E.) 统计 量 (Waldx”) 自由 度 (4/) 显著 性 水 平 (P) 发 生 上 比率 
Independent variable Regression Stand error Statistics Degree of Significant (OR) 

P coefficient freedom level Odds ratio 
高 程 Elevation —0.003 <0.001 900.270 1 <0.00 0.997 
坡度 Slope —0.006 0.002 5.545 1 0.019 0.994 
照射 坡 向 Trasp -0.640 0.033 367.300 1 <0.00 0.527 
十 计 E 
坡 位 指数 —0.008 <0.001 462.390 1 <0.00 0.992 
Slope position index 
地 形 起 伏 度 Relief -0.029 0.001 2 045.482 1 <0.00 0.972 
年 平均 气温 
Annual average 0.220 0.026 70.173 1 <0.00 1.246 
temperature 
寒冷 指数 Coldness index 0.140 0.013 121.506 1 <0.00 1.150 
干燥 度 Drought index -1.192 0.305 15.262 1 <0.00 0.304 
常数 Constant 7.590 0.919 68.185 1 <0.00 1 978.469 


HL (homsmer-lemeshow) = 2.532. P=0.98. 
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3 ” 井 队 县 影响 耕 资源 的 生态 要 素 栅 格 图 
Fig.3 Grid maps of the key ecological factors impacting cultivated land resources of Jingxing County 
A: 高 程 (m); B: 坡度 (*); C: 照射 坡 向 (无 量 纲 ); D: 地 形 起 伏 度 (m); E: 坡 位 指数 (m); F: 年 平均 气温 (C); G: 年 平均 降水 (mm); H: 
年 活动 积温 ('C); I: 年 生物 积温 ('C); J: 温暖 指数 (C); K: 寒冷 指数 (CC); L: 湿润 指 数 (mm C7); M: 干燥 度 (无 量 纲 )。A: elevation (m); 
B: slope (°); C: trasp (non-dimensional); D: relief (m); E: topographic position index (m); F: annual average temperature (‘C); G: annual 
average precipitation (mm); H: active accumulate temperature (‘C); I: biological cumulative temperature (C); J: warmth index (C); K: 
coldness index (C); L: humidity index (mm:°C !); M: drought index (non-dimensional). 
影响 山区 耕地 资源 潜力 的 8 个 关键 生态 要 素 的 贡献 
率 排序 为 : 地 形 起 伏 度 > 高 程 > 坡 位 指数 > 转换 坡 向 > 
寒冷 指数 > 年 平均 气温 > 干燥 度 > 坡度 ， 所 选 的 生态 
要 素 中 地 形 要 素 的 贡献 率 大 于 气象 及 气象 相关 要 素 


的 。 因 此 , 评估 耕地 太行 山 丘 陵 区 耕地 资源 潜力 时 ， 


地 形 要 素 配 置 的 分 析 较 气象 及 气象 相关 要 素 的 配置 人 
更 为 重要 。 
由 模型 参数 中 的 回归 系数 5( 表 3) 分 析 得 出 各 
个 生态 要 素 与 山区 耕地 资源 潜力 的 正 负 相 关 关 系 。 Ger 
8 个 影响 山区 耕地 资源 的 关键 生态 要 素 中 ， 除了 年 eg 
平均 气 气温 昌 和 寒 准 全 指数 ， 其 余 变量 因子 都 是 负 的 Cultivated land potential 
回归 系数 ,说 明 随 着 年 平均 气温 和 寒冷 指数 的 增 0.99 
加 ,山区 耕地 资源 潜力 呈现 增加 的 趋势 ,而 其 余 自 一 
变量 因子 的 增加 会 引起 山区 耕地 资源 潜力 呈现 减 0.00 
小 的 趋势 。 
变量 因子 的 改变 引起 山区 耕地 资源 潜力 变化 图 4 基于 二 分 类 Logistic 回归 模型 的 井 取 县 耕地 资源 潜 
可 量化 , 由 表 3 发 生 比 率 (OR) 可 知 ， 除 了 干燥 度 改 | 。” ” ”” 力 分 布 图 


Fig.4 Distribution map of cultivated land potential of 
变 1 个 单位 量 ,山区 耕地 资源 潜力 变化 0.3 倍 ， 其 余 Jingxing County based on the binary Logistic regression model 
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变量 因子 改变 1 个 单位 量 ,山区 耕地 资源 潜力 变 
化 1 倍 左右 。 
2.2 ”耕地 资源 潜力 估算 

将 表 3 中 的 回归 系数 带 入 到 方程 (1)， 得 到 
Logistic 回归 模型 ， 根 据 模型 利用 GIS 运算 生态 要 素 
栅 格 图 ， 得 到 井 取 县 耕地 资源 潜力 分 布 图 (图 4)。 由 
图 4 可 以 看 出 , 井 取 县 耕地 资源 潜力 在 0.5 以 下 地 区 
所 占 面积 极 少 。 对 于 该 模型 ， 人们 一 般 选 择 的 划分 
阅 值 为 0.5, 但 是 从 生态 安全 的 角度 ， 山 区 土地 资源 
成 为 耕地 资源 的 条 件 比较 严格 0 1， 因 此 选择 0.9 作 
为 阀 值 来 划分 研究 区 域 能 够 成 为 耕地 的 地 区 。 耕 地 
资源 潜力 大 于 0.9 的 分 布 区 如 图 5 所 示 。 分 析 可 得 ， 
基于 二 分 类 Logistic 回归 模型 ， 井 降 山区 土地 资源 
能 够 成 为 耕地 的 面积 为 60 400 hm*, 应 用 SPOT5 影 
像 分 析 得 到 井 降 县 的 耕地 面积 45 600 hm ， 故 井 降 
县 还 具有 14 800 hm 的 耕地 后 备 资源 ,， 相当 于 现 有 
耕地 面积 的 32.5%。 由 此 可 知 ， 井 取 县 上 县 有 较 多 的 耕 
地 后 备 资源 。 


国 国 现存 耕地 Existing cultivated land 
非 耕 地 Non cultivated land ,Vu 
潜在 耕地 Potential cultivated land ” ee 


5 基于 二 分 类 Logistic 回归 模型 井 降 县 后 备耕 地 资源 
分 布 
Fig.5 Distribution map of cultivated land reserve of Jingxing 
County based on the binary Logistic regression model 


3 结论 与 讨论 

本 文 从 生态 角度 出 发 ,选择 13 个 山区 耕地 资源 
潜力 的 生态 要 素 , 利用 二 分 类 Logistic 回归 模型 分 
析 ， 得 到 8 个 关键 生态 要 素 , 分 别 是 地 形 起 伏 度 、 高 
程 、 坡 位 指数 、 转 换 坡 向 、 骞 冷 指数 、 年 平均 气温 、 
干燥 度 、 坡 度 。 基 于 建立 的 Logistic 回归 模型 ， 利 用 
Arcgis 软件 得 到 井 降 县 耕地 资源 潜力 分 布 图 ， 经 分 
析 得 出 井 取 县 还 上 县 有 14 800 hm 的 后 备耕 地 资源 ， 
相当 于 现 有 耕地 面积 的 32.5%, 为 以 后 井 取 县 耕地 
后 备 资源 的 开发 利用 提供 了 一 定 的 理论 依据 。 


在 山区 ,高程 和 坡度 是 影响 耕地 资源 分 布 的 最 

影响 因素 ”1 时。 在 丘陵 区 地 形 起 伏 度 和 高 程 是 
影响 丘陵 区 耕地 资源 潜力 的 两 大 最 主要 的 生态 要 素 ， 
其 贡献 率 大 于 70%。 因 此 ， 高 程 是 影响 丘陵 区 耕地 
资源 分 布 的 主要 因素 之 一 , 这 与 前 人 研究 结果 一 
致 。 但 是 ， 丘 陵 区 坡度 并 不 是 影响 丘陵 区 耕地 资源 
分 布 的 最 主要 的 因素 , 而 地 形 起 伏 度 则 是 影响 耕地 
资源 分 布 的 主要 因素 之 一 。 其 原因 是 太行 山 丘 陵 区 
坡度 较 缓 , 80% 的 地 区 坡度 <25°( 图 3B)。 

张 迪 等 "指出 耕地 后 备 资源 中 ,面积 最 大 的 是 
荒草 地 ， 这 部 分 荒草 地 大 多 分 布 在 山地 丘陵 区 。 张 
贵 军 等 (1 按 不 同 地 貌 类 型 将 河北 省 划分 为 8 个 类 型 
区 ， 认 为 太行 山 丘 陵 区 耕地 后 备 资源 潜力 较 大 。 本 
研究 结果 进一步 证 明了 前 人 的 研究 结论 。 但 是 ， 本 
研究 所 选择 的 决定 耕地 资源 分 布 的 生态 要 素 中 没有 
包括 土壤 因子 , 而 土壤 是 决定 一 个 地 区 是 否 能 够 成 
为 耕地 的 决定 性 条 件 0 ”4 。 由 于 山区 土壤 数据 较 
难 获 得 ， 且 缺少 研究 区 土壤 理化 性 质 以 及 土 层 厚度 
的 空间 化 数据 。 另 外 , 根据 我 国 1 : 400 万 的 土壤 类 
型 图 ， 井 降 县 土壤 类 型 较 一 致 ， 全 部 为 淋 溶 褐 土 。 因 
此 ， 本 文 建立 井 阴 县 耕地 模型 中 , 没有 考虑 土壤 因 
素 。 对 于 更 大 区 域 的 研究 中 , 为 了 模型 的 准确 性 ， 应 
该 加 入 土壤 因素 。 
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